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We present a detailed analysis of the radio synchrotron emission induced by WIMP dark matter 
annihilations and decays in extragalactic halos. We compute intensity, angular correlation, and 
source counts and discuss the impact on the expected signals of dark matter clustering, as well as 
of other astrophysical uncertainties as magnetic fields and spatial diffusion. Bounds on dark matter 
^— ^ ■ microscopic properties are then derived, and, depending on the specific set of assumptions, they are 

' competitive with constraints from other indirect dark matter searches. At GHz frequencies, dark 

matter sources can become a significant fraction of the total number of sources with brightness 
^ ■ below the micro Jansky level. We show that, at this level of fluxes (which are within the reach of 

the next-generation radio surveys), properties of the faint edge of differential source counts, as well 
as angular correlation data, can become an important probe for WIMPs. 
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PACS numbers: 95.35.+d,95.30.Cq,95.85.Bh 



I. INTRODUCTION 



The nature of non-baryonic Dark Matter (DM) is one of the most puzzling mysteries of modern Cosmology (see 
e -g- EL f° r reviews). A non-gravitational detection of dark matter would be a fundamental step toward the 
comprehension of its properties. This could be achieved in several ways, the most promising and popular methods 
C$ ■ are direct searches in underground experiments, DM searches at colliders, and indirect detection methods. The latter 
techniques are based on the detection of the products of DM annihilations or decays, notably neutrinos, antimatter 
■ and photons. Indirect searches are particularly promising if dark matter is in the form of Weakly Interacting Massive 
1 Particles (WIMPs). In the early Universe, these particles decouple non-relativistically from the thermal plasma and 
r**"- . inherit the correct relic abundance for a thermally averaged annihilation cross section close to (av) — 3 x 10 -26 cm 3 s~ 1 . 

Annihilations of WIMPs inside galaxies are able to inject large amounts of relativistic electrons and positrons, which in 
*n | turns interact with the galactic magnetic fields producing synchrotron radiation. This is a generic prediction of WIMP 
NJ . • models, except for the peculiar case of WIMP candidates annihilating only into neutrinos. For galactic magnetic fields 
| of the order of /iG and electrons energies below about 10 GeV, the synchrotron emission falls at frequencies around 
and below the GHz , i.e. in the radio band. Therefore, radio observations arc appropriate tools for indirect WIMP 
searches and they have been extensively discussed in the context of WIMP annihilations inside our galaxy jHT|. 
In particular, low frequencies (< 1 GHz) data severely constraint WIMP masses < 10 GeV for leptonic annihilation 
modes [l5l |. 

Here we focus instead on the radio emission produced by DM annihilations or decays in extragalactic halos. A 
previous analysis along this direction can be found in Ref. [16[ . Recently, we have shown that for realistic assumptions 
d \ on the DM clustering and magnetic fields and for a "thermal" annihilation cross-section (i.e. (av) = 3x 10~ 26 cm 3 s _1 ), 
the DM signal could account for a significant fraction of the extragalactic-radio emission inferred by the ARCADE 2 



Collaboration [17J . This constitutes a further motivation to seriously consider DM searches with extragalactic radio 
observations. 

In this work, we extend the studies in [l6l Il7j by extensively analyzing the impact of astrophysical and particle 
physics uncertainties on the DM extragalactic radio emission which is then compared to observational data. We 
focus on three observables: intensity, differential number counts of sources and angular correlations. By comparing 
predictions with current data, constraints on the DM annihilation/decay rate and mass are derived. We also compare 
the DM signals with (data-driven) expectations for astrophysical sources, such as radio loud active galactic nuclei 
and star forming galaxies, and show that source counts and angular correlation data could be particularly relevant to 
distinguish the DM contribution from the astrophysical ones. 
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It is interesting to note that cosmological radio emission induced by WIMP DM decreases very rapidly with redshift 
(so the main contributions come from structures at z -C 1), and that DM halos can substantially contribute to the 
counts of sources at brightness below the /xJy level. Therefore, a specific property of the extragalactic WIMP-induced 
source population is to peak at low redshift and low brightness. This hypothesis can be tested by telescopes with flux 
sensitivity at /xJy level, which will be reached by upcoming radio experiments, in particular by the Square Kilometer 
Array (SKA). 

The paper is organized as follows. In Section III Al we review the present status of radio observations, focusing 
on those particularly relevant for DM searches. We discuss data on intensity, angular correlations and counts of 
sources. The formalism needed to compute the DM and astrophysical signals is presented in Section III Bl Section 
IIII Al shows our results for both DM and astrophysical models. Our estimates are first derived focusing on few specific 
DM benchmark cases while a detailed analysis about the impact of astrophysical and particle physics uncertainties is 
presented in Section IIII Bl The bounds on the microscopic properties of DM are computed in Section IIII CI and we 
derive our conclusions in Section [IVI 



Although the history of radio observations is quite long and rich, only few maps have sufficient angular resolution 
and sky coverage to allow estimates of the isotropic extragalactic background. In particular, they need to be suitable 
to reliably separate the large-scale (astrophysical) Galactic component, which is a delicate procedure. 

Recently, the ARCADE 2 (Absolute Radiometer for Cosmology, Astrophysics and Diffuse Emission) pro ject (l8j 
reported the measurement of the extra-galactic sky temperature at frequencies ranging from 3 to 90 GHz |l9l . l20j . At 
high frequency, the radiation is largely dominated by the CMB, so the estimate of other extra-galactic components 
becomes unfeasible. At frequencies down to GHz the CMB is still dominant but its contribution can be easily 
subtracted since it follows a black-body (BB) spectrum and we know its BB temperature with per-mil precision. At 
lower frequencies, extra-galactic sources dominate. 

In the ARCADE data, the diffuse Galactic foreground is separated using two independent estimators (namely, a 
co-secant dependence on Galactic latitude and the correlation between radio and atomic line emissions), which agree 
well between each other [2l|. In Ref. similar analyses have been performed on high-latitude patches in past 
surveys at 22, 45, 408, and 1420 MHz. We consider the estimates derived with this procedure as a conservative 
constraint for the extra-galactic signal since they might include a contribution from a (isotropic) Galactic emission 
uncorrelated with th gas (e.g., from the DM halo). Data-points shown throughout the paper are taken from Table 4 
in Ref. [l9| by converting the thermodynamic temperature To to the brightness temperature Tb by means of the 
usual procedure: Tb = To • x/(e x — 1) with x = hv/(kTo). 

The expected signal of astrophysical extragalactic sources derived from models which fit source counts (see next 
Section) lies well below the data: this points towards the existence of a previously disregarded component. The 
extragalactic (or Galactic) DM emission can account for this excess and this possibility has been explored in Ref. [ljj ■ 

We will also consider estimates of the cosmic infrared background derived in Ref. [22| (with CMB subtracted) and 
based on FIRAS and DIRBE datasets. 



The angular distribution of sources in the sky is a powerful probe of large-scale clustering. Probing the latter 
(which includes disentangling contributions from different redshift slices) at radio frequencies is rather hard, but, 
nevertheless, wide-area radio surveys are a useful tool, since they can test clustering on the largest scales. Recently, 
observational campaigns (e.g., FIRST [23| and NVSS [24||), have been exploited for this purpose. 

Angular distributions are often reported in terms of the two-point angular correlation function lo{6) or the an- 
gular power spectrum Cg. Theoretically, such two estimators are fully equivalent and simply related by a Legendre 
transformation: 



II. EXTRA-GALACTIC RADIO SIGNALS 



A. 



Observational data 



1. Total intensity of cosmic background 



2. Angular correlation 




(1) 
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where N is the number of sources, AQ is the survey area, and Pg, are the Legendrc polynomials. On the other hand, 
Ci and u>(0) quantify different properties of the angular distribution and so they are not observationally equivalent 
(i.e., to perform the transformation we would need data at all angles, which is of course unrealistic). Naively, Ce is a 
better probe at large angles (provided a sufficiently large sky coverage), while w, which suffers from Poisson noise on 
large scales, is a superior tool for small angular separations. 

For what concerns the angular two-point correlation function of extragalactic sources, the state of the art bears 
on flux thresholds of few mJy and angular separation of (fraction of) arcmin. Since the DM contribution is faint, 
we aim at analyzing correlation functions with flux threshold as low as possible. We consider data from the FIRST 
survey [25[ at 1.4 GHz analyzed in Ref. [26| (their Fig. 3), which concerns clustering of sources with S > 3 mJy. For 
our purposes, although other measurements could be considered (including surveys of VLSS at 74 MHz [271 ]. WENSS 
at 325 MHz (H, SUMSS at 843 MHz [H], and NVSS at 1.4 GHz J^), they do not add further insights since they have 
similar (or worse) flux thresholds. We anticipate already here that the number of DM sources with mJy brightness at 
GHz frequencies is very low for typical WIMP models, thus the clustering tested by those surveys is not really relevant 
to DM. To partially overcome such issue one can consider measurements of angular correlation functions in small but 
deep fields (e.g., at 1.4 GHz in [3(3] with S > 0.2 mJy). However, in this case the statistics is much poorer and this 
makes any interpretation significantly harder. As we will discuss later, flux thresholds at /xJy level would be in order 
to possibly probe realistic WIMP scenarios, and this can be actually achieved only by future surveys with SKA. The 
angular power spectrum CV of the radio source distribution can be inferred (under some assumptions and using Eq. [T]) 
from the angular correlation function. This has been performed in Ref. [311 ] for multipoles below £ ~ 100 analyzing 
NVSS data at 1.4 GHz. However, large statistical errors and, again, high flux thresholds limit their usefulness for 
constraining WIMP models. 

Radio C( of anisotropics are mainly extracted from sky-maps of single-dish radio telescopes at low-frequency and 
from CMB-orientcd experiments at high-frequency. The firsts have typically low angular resolution, and provide 
estimates of Ci up to £ of few hundreds, while the Galactic foreground dominates up to I ~ 100. We consider ang ular 
power-spectra at different frequencies derived in Ref. 132] with estimates of source contribution as in Ref. [U [Hj]. 
We derive a band which brackets uncertainties related to the disentanglement between extragalactic and galactic 
foreground anisotropics. 

At CMB frequencies, the CMB fluctuations dominates over extra-galactic sources except for v > 100 GHz and 
£ > 3000 (where unfortunately IR sources constitute the main contribution). Wc consider recent data at 150 GHz 
from SPT [35| (similarly, the analysis of ACT data [3|| leads to analogous conclusions). 

Therefore, at present, measurements of radio angular correlation do not provide very powerful constraints on non- 
gravitational DM signals. On the other hand, our analysis is mainly oriented towards understanding anisotropics 
signatures of WIMPs and how they could be tested by future surveys. 

3. Source number counts 

Source number counts arc an important tool to understand properties of different radio populations. Although 
counting sources of a given brightness may sound rather trivial, it is well known that obtaining a population count 
from a sky survey is not an easy task, and different telescopes (endowed with, e.g., different angular resolution and 
sensitivity) introduce different biases that need to be properly accounted for. A wide consensus has been reached 
nowadays on the interpretation that radio counts are considered to be dominated by radio-loud AGNs down to mJy- 
levels, while star-for ming galaxies (and possibly radio-quiet AGNs) take over at fainter flux-density. For a recent 
review, see, e.g., Ref. [371 ] . 

The number of sources per unit area N is often reported in terms of the differential counts dN/dS (where S is the 
flux density) and normalized to the Euclidean case, i.e., multiplied by S 5 ^' 2 (a uniform source distribution in a static 
Euclidean universe leads to N oc 5 -3 / 2 , since the number of sources is proportional to the volume, while the flux 
density scales as the inverse distance squared). 

A collection of differential source counts ranging from 150 MHz to 20 GHz is provided in Ref. [37j (see references 
therein for each dataset), and will be exploited in this work. 
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B. Theoretical models 



Considering an extragalactic source population with number density n and with luminosity C depending on a certain 
quantity Q, the total isotropic intensity per solid angle at a given frequency v is given by [la. l38j: 

d2y [in dn in d^ + z)C{E = E v {\ + z) lZ ,Q)e~ T ^ 



4?r d L (z) 



2 



where V is the comoving volume, f2 is the solid angle, z is the rcdshift, E is the emission energy, dj, is the luminosity 
distance, H is the Hubble rate, and r is the optical depth. Q will be M (halo mass) for DM and £„ (luminosity at 
a reference frequency vq) for astrophysical sources. 

The differential number counts of sources per solid angle can be written as: 

dN f e-^UUz) 4 dn,„ dQ 

is; = 47TC J dz (i + zrH(z) dQ iQ > z) dc;> (3) 

where S v = (1 + z) £„, /(4it d L (z) 2 ) and v' = (1 + z) u. 
Eqs. ^ and ([3]) can be simply linked through: 

s ° dN 

ds k s ' < 4) 

where So is a minimum (telescope-dependent) apparent source-brightness above which sources can be detected by 
the telescope and so subtracted from the total isotropic flux (we will mostly refer to l v as the contribution from 
unresolved extragalactic sources). 

The angular power spectrum of cosmic background anisotropies is fl6l . [38| : 

C* £ = C} h + Cj h (5) 

r f p - 2t ( z ) r dn 

cr = ^/^^/^^ (6) 

n 2 

, (7) 



r f p 2t ( z ) r r]n 

C « = JA^Y J dz d L {zYH{z) Plm{Kz) J W-jQ(Q,z)£(E = E v (l + z),z,Q)b(Q,z)\u(k\Q)\ 

where u is the Fourier transform of the spatial emission density profile of the source (with the latter normalized such 
that its volume integral is equal to unity) , and b is the linear bias factor that relates the source clustering to the matter 
clustering. Label "1ft," denotes correlations between particles within the same halo, while "2/i" describes correlations 
between particles in two distinct halos (and since wc are interested in the two-point correlation, there are no further 
terms). 

The Limber equation provides the relation between number counts and angular power spectrum or correlation 
function of a given class of objects [H, E3|- ^ n particular, we will make use of (see, e.g., Ref. |4lj): 

f°° „ , f°° , „ , Jo(kOr) { 1 dN\ 2 dz 

where P(k, z) is the 3D power-spectrum, Jo is the zero-order Bessel function, dz/dr = H(z)/c in a flat universe, and 
dN/dz can be derived from Eq. ([3]). Note also that, in the idealized case of uniformly distributed point-sources, Ci 

does not vary with i (i.e., it is a Poisson noise term) and is simply given by Ce = J Q ° dS dN/dS S 2 . 

Absorption along the line of sight (described through the optical depth r) is negligible at radio frequencies above 
~ 10 MHz and will be disregarded. Possible absorption effects within the source (e.g., synchrotron self-absorption) 
are instead included in £; however, for the (non-compact) DM sources considered here, they are negligible as well. 



Dark matter source 



We consider the DM to be in form of particles (with mass M x ) which can annihilate (with annihilation rate {cr a v)) 
or decay (with decay rate t^ 1 ) into (among other products) electrons and positrons. The latter in turn emit photons 
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through synchrotron radiation. Up to the injection of e~/e + and for a given DM model, the process is completely 
set by the amount of DM present in the source, which implies Q = M, where M is the DM halo mass. On the other 
hand, some of the quantities involved in the diffusion and energy loss of e~/e + after production (e.g., magnetic field) 
could be not strictly correlated with the mass of the source, but rather to properties of host galaxies, merger history, 
etc. As a first approximation, however, their possible variations from source to source can be reasonably associated 
only to halo mass and redshift, and so, for simplicity, we will assume Q = M . The luminosity C of DM sources is 
given by 



C h a h (E, z, M) = E ^ [(1 - /) p(M, r, z)] 2 for the host halo (9) 

C*{E,z,M) = E Y^j™ B dM s ^(M s J,M) J*" d 3 r 3 ^p 2 s (M 8 ,r s ,z) forsubhalos (10) 



in the annihilating DM case, and by 

C d (E, z, M) = —— / P {M, r, z) (11) 

Td m x j at, 

in the decaying DM case. We introduced p which is the total DM density profile (cut at virial radius R v ), and p s 
which is the density profile of a single clump. Note that in Eq. [10]thcrc is no dependence on the spatial distribution of 
clumps within the host halo. This is not relevant when calculating the total isotropic intensity, while the information 
is encoded in u for the angular power spectrum. The spatial profile of the main halo is given by pu = (1 — f) p. 
Thus / provides the fraction of total halo mass given by substructures with mass function dn s /dM s (normalized such 
that / dM s dn s /dM B M s = f M). The minimum subhalo mass is M* ut , and dNi/dE = dNi/ dE{E,r, z). We aim at 
providing equations suitable for multi-wavelength analyses and the subscript i denotes different emission mechanisms; 
more precisely 



dN~, 


dN y 


dE 


dE ' 


dN svn jc 




dE 


= 2 y 



(12) 



dE' ^v™* 10 . n e ) f or radiative emission . (13) 
E 

where h e (r,E) = n e /A, with A a = (av)/2 ■ (p/M x ) 2 and Ad = t^ 1 p/M x , and n e being the electron/positron equilib- 
rium number density obtained solving a transport equation for e~/e + injected by DM with an energy spectrum set 
by dN e /dE (see, e.g., review in Ref. [42J). The function dNj/dE with j = 7, e describes spectra of annihilation/decay 
into gamma-rays and electrons/positrons, respectively. The radiative emission power in the synchrotron case is given 
by: 

P syn {r, E, v) = ^—j B{r)F{v/ Vc ) (14) 
ui e c 

where m e is the electron mass, the critical synchrotron frequency is defined as v c = 3/(4 7r) • ce/(m e c 2 ) 3 B(r)E 2 , and 
F(t) =t J. dzK 5 / 3 (z) is the function setting the spectral behavior of synchrotron radiation. For the inverse-Compton 
(IC) process we have: 



Pic(r,E,u) = chv J de n 7 (e, r) er(e, u, E) (15) 

where e is the energy of the target photons, n 7 is their energy spectrum, and a is the Klein-Nishina cross section. 

The quantity u in Eqs. |6] and [7] is the Fourier transform of the spatial emission density profile U of the source 
normalized to unity (i.e., J d 3 r U = 1). In the annihilating case, for the main halo, U is proportional to p^, and turns 
out to be simply given by A a ■ E/v ■ dNi/dE up to the normalization factor. In the substructure case, we consider a 
scenario in which U follows from p^ of the parent halo (biased) , and an antibiased case modeled from results in [43[ 
(which are assumed to hold for any M ) where the subhalos are more preferentially located in the outer part of the 
parent halo. With this treatment, clumps are effectively considered as point sources with spatial probability set by U, 
thus disregarding the actual dumpiness of subhalo distribution which however would affect only angular correlations 
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at extremely small angles. 1 

For decaying models, the spatial emission density profile of the source is given by ■ E/v ■ dNi j dE with p = Ph + p s 
in the definition of Ad- 

Neglecting spatial dependencies in dNi/dE, this term goes out of the volume integral in Eqs. [5J EH and [TO This 
implies also that n e is spatially constant which can be satisfied only if in the transport equation for n e there is no 
spatially dependent term (except for the source term itself which in this case encodes all the spatial information). 
In particular, under this approximation, electrons and positrons radiate at the same place where they are injected, 
without undergoing spatial diffusion. In this scenario the solution of the transport equation can be written as: 

n e (E, z) = dE' A^(E') , or equivalent^, n e {E) = dE' ^{E') . (16) 

where ti is the time-scale associated to energy losses. We consider this scenario as our first approximation in Sec. MI Al 
(while impacts of possible spatial dependencies of the magnetic field and spatial diffusion of e + /e~ in the source will 
be discussed in Sees. MI B 31 and MI B 4ft . 

With these assumptions and disregarding the contribution of substructures, Eq. ([2]) reduces to (see Ref. [Hj]): 

c.E 2 (a a v) fn xPc \ 2 f , (l + z) 3 A 2 (z) dNi t(z ) r • ™ r 

"'- = i^VUih) J* - 4) -d " for—tog DM, (17) 

- ^-MI dz wm^" f «*»™ DM > < 18 > 

with p c being the critical cosmological density. The quantity A 2 (z) = (<5 2 ) (with d = p/pcDM where pcdm is the 
mean DM density and we will consider pcdm = Px — Pc) represents the clumping factor, and encodes all the 
effects due to the DM clustering. It is extensively described in, e.g., Ref [45[. Effects related to the DM halo mass 
function, the density profile within each halo, and the minimum halo mass in the clumping factor, are discussed in 
Sec. MI B II 

Plugging Eqs. (|10j) and (fTTj) into Eqs. (JB]) and (JT]), the angular power spectrum of DM can be derived. Under the 
same assumptions that lead to Eqs. (jTTJ) and (|18[) . one can proceed to analogous simplifications, and again equations 
found in the prompt gamma-ray emission case (e.g., in Ref. [HI, H^|) can be applied to the radio case by simply 
replacing dN^/dE with dN syn /dE. More precisely, in such a simplified scenario, one finds: 

Q=(^)"(^) / 777^ 1(1 + *>'§)' = </'.«» for— ingDM, (19) 

°> = [i^fk ) J hww p,Ak = " r - =) c " n " for dcc * ymg DM ■ (20) 

where Pf t is the 3D power spectrum of the Fourier transform of /j, where f a = S 2 — (S 2 ) in the annihilating case, 
and fd = d— (5) in the decaying case (thus Pf d is equal to the 'standard' non-linear matter power spectrum; for the 
latter, see, e.g., the review in Ref. [13|). The expression of Pf a can be easily derived adding Eqs. [6] and [7] (for further 
details, see Ref. [38[) which leads to: 




P fa {k,z) = P} } : + P^= J dM-^\u(k\M)\ 2 +y dM^b(M,z)\u(k\M)\ 
with u(k\M) being the Fourier transform of p 2 {r\M)/p 2 . 



2 

Hit 



P lm (k,z), (21) 



2. Astrophysical sources 

The radio cosmic background from astrophysical sources is thought to mainly come from two distinct populations, 
radio loud active galactic nuclei (AGN) and star-forming galaxies (SFG) (for a recent review on radio sources, see, 



1 Our approach is fully equivalent to the description of Ref. |44j . by rewriting their Halo Occupation Distribution as (N\M) = 
/ dM s dn/dM s f d 3 r s p 2 s which implies C s a h = E z (a a v)/(2M 2 ) dN t /dE (N\M). 
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FIG. 1: Total intensity. Extragalactic radio (left panel) and gamma-ray (central panel) backgrounds for benchmark Bl (blue) 
and B2 (orange) WIMPs in the annihilating (solid) and decaying (dotted) cases (see text and Table[T]for details). For the radio 
case, main astrophysical source contributions (black-solid), CMB (black-dotted), and the best-fit of ARCADE data (black- 
dashed) are also shown [53]. The right panel shows contributions at different redshifts (i.e., f* dz' dl /dz') for emissions 
given by synchrotron radiation at 1 GHz, inverse Compton on CMB at 1 MeV, and 7r°-decay at 1 GeV, in the benchmark case 
Bl (all normalized to unity). 

e.g., Ref. [13] )■ The first dominates source counts at high fluxes, while the latter takes over below ~mJy level. Below 
mJy level, radio quiet AGNs can also provide an important contribution (4^|, which is nevertheless subdominant, and, 
for simplicity, we will disregard it. The luminosity of astrophysical sources is typically written in terms of luminosity 
at a reference frequency vq where the most robust luminosity function models can be built (i.e., at frequency where 
more data are present) and including a running to the frequency of interest. We will adopt this method and take 
Q = £„ with Vq = 1.4 GHz and Vq = 60 fim for AGN and SFG, respectively. 

The comoving luminosity function <j)(C,z) = dn/dlogC of AGN at 1.4 GHz is taken from Ref. [49j| where cf> = 
no/[(£/L*) a + (£/L*) b ] (for values and definition of no, a, b, and L* parameters we refer to their Table 1). The 
luminosity at other frequencies is obtained assuming a spectrum £ oc v~ a . The model includes two flat-spectrum 
populations (a = 0.1) with different evolutionary properties (BL Lacs and Flat Spectrum Radio Quasars) and a 
steep-spectrum population (a = 0.8). 

Star forming galaxies are typically described by means of their properties in the infrared (IR) band. Although the 
IR emission may not be generated by the same physical mechanism of radio emission, radio sources are associated 
to young stars and so tightly correlated to the star formation rate probed at IR frequency. Empirically, the far 
IR-radio correlation is, at least for high-brightness and low redshift, well established, and follows a nearly linear 
relation [37j (that we will assume to be redshift independent): Ciaghz = 1-16 • 10~ 2 Lb / '[Lb /Ceo^m + (Lbf^-sopm) 3 ' 1 ] 
with Lb = 8.8 • 10 29 ergs _1 Hz _1 . The luminosity at frequencies other than 1.4 GHz is again obtained through 
£ oc v~ a with a = 0.75 (from typical spectral index of optically thin synchrotron sources). We consider two SFG 
populations which are expected to be reservoirs of cosmic-rays and so bright in synchrotron emission, starburst 
galaxies and normal late-type spiral galaxies. Their luminosity function dn/dC has been modeled following Ref. [oOj 
with parameters derived from, respectively, 'warm' and 'cool' IRAS galaxies in Ref. [5lj ]. Our modeling of AGN and 
SFG agrees well with number counts and angular correlation data. 

To compute the angular power spectrum of anisotropics generated by astrophysical sources, we would need to 
model the spatial profile of luminosity within the source (to get its Fourier transform u in Eqs. ([6]) and ([7])). This 
is not as simple as in the DM case (where it follows a universal profile) and in the following we treat astrophysical 
sources as point sources (u = 1). Although this is a crude approximation, it does not affect results up to t ;§> 10 3 
or 9 <C 0.1° (see, e.g., Ref. [29[, where it is found that multiple components of radio galaxies become important in the 
correlation function only below few arcmin scales). With this assumption the form of each of the two terms of the 
angular power spectrum in Eq. ([5]) is known, with C\ h ~ const, and the ^—dependence of Cf h set by P hn . However, 
their absolute size depends on £ and on luminosity integral extrema, so their ratio and the form of total Ci is model 
(and telescope) dependent. 
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Name 


Mass 
[GeV] 


(a a v) [cm ,! s 1 ] 
annihilating case 


r [s] 
decaying case 


Dominant 
final state 


Bl 


100 


3 ■ 10 -26 


4 ■ 10 28 


b-b 


B2 


10 


3- 10- 26 


5 ■ 10 27 


M + - M~ 



TABLE I: Benchmark DM models. 



III. RESULTS 
A. Benchmark cases 

In order to have a feeling about the expected DM signals and how they compare to astrophysical emissions, we first 
show results for few benchmark cases. We consider a neutralino-likc DM candidate (model Bl) with mass M x = 100 
GeV, 'thermal' annihilation rate (a a v) = 3- 10 -26 cm 3 s _1 and final state of annihilation into quarks b—b (and similarly 
a decaying case with r = 4 • 10 28 s). Motivated by recent claims from direct detection experiments and by a possible 
DM interpretation of the ARCADE excess [13], we introduce also a DM candidate (model B2) with mass M x = 10 
GeV, again with 'thermal' annihilation rate, but annihilating into the leptonic final state [i + — fj,~ (to enhance the 
synchrotron signal), and an analogous decaying case with r = 5 • 10 27 s. 

We adopt the halo mass function from Ref. [52J , the concentration of halos from Ref. [53[ , and the DM distribution 
in halos following a NFW profile (lij . No substructures are considered in this benchmark scenario and the minimum 
halo mass is set to M cut — 10 5 M Q . The magnetic field is assumed to be constant in space and time with B = 10 /iG, 
and e + /e~ are assumed to radiate at the same place where they are injected. Each of these assumptions will be 
discussed and relaxed in the next Sections. 

Throughout the paper, we will stick to a 6-parameter ACDM model derived by the WMAP team combining CMB 
data with Baryon-Acoustic-Oscillation (BAO) and Hq measurements [55j] 2 . Parameters relevant for our purposes are 
fi m = 0.27, il A = 0.73, n CD M = 0.23, h = 0.70, n s = 0.96, and cr 8 = 0.8. 

In Fig. QJ,, we plot the total intensities. Note that the DM contribution in those benchmark scenarios is roughly at 
the same level of astrophysical emissions. The extra-galactic sky temperature reported by the ARCADE 2 collabora- 
tion [l9| cannot be explained by ordinary models of AGN + SFG. In Ref. [13], we discussed the possibility of fitting 
ARCADE data in terms of a DM-induced synchrotron emission. WIMP models can actually account for the excess, 
with viable scenarios being only slightly more optimistic than the benchmark cases considered here [13]. Comparing 
Fig. [TJi with Fig. [T}), one can see that radio extragalactic data are more constraining than the 7-ray counterpart for 
WIMP models with dominant annihilation/decay channel into leptons (by roughly one order of magnitude in the B2 
case); on the contrary, for WIMP models with significant 7r° production (as in the case of dominant hadronic final 
states Bl), gamma and radio data have roughly the same constraining power. 

We have been considering energy losses due to synchrotron radiation and to IC scattering on CMB only. The 
effect of the electron/positron interactions with the interstellar radiation field (ISRF) is difficult to model within this 
formalism since it is mostly related to the properties of the astrophysical sources hosted within each halo, rather 
than to the halo mass. However, the IC loss associated to the ISRF is expected to be generically subdominant in 
the computation of extragalactic signals. To understand this point, we included energy losses on a simple ISRF with 
constant density of lev cm" 3 in all halos at any M and z. This density is roughly the mean optical ISRF density 
of the Milky Way and can be considered as an overestimate of the actual field since it is known that either lager 
(e.g., cluster) and smaller (e.g., dwarf spheroidal galaxies) structures have much lower density. The effect on the total 
intensity curves in Fig. [TJl is an overall depletion of ~ 15%. We therefore assume energy losses due to IC on ISRF to 
be subdominant, and we disregard them in the rest of the paper. 

In Fig. QJ,, the redshift distribution of main radio, X-, and 7-ray signals is shown. In the case of inverse Compton 
scattering on CMB, it is completely set by the evolution of the DM spatial distribution. Indeed, the dilution of CMB 
energy density due to the expansion of the Universe affects both energy losses and IC power, having basically no net 
effects on the emission. There can be a further suppression in the emission at low z if synchrotron losses take over at 
late times. Moreover, the location of the IC peak is roughly redshift independent (and so the energy of the emitted 
electrons), since the factor 1/(1 + z) arising from emission to the observer is exactly compensated by the increase 



2 Note, however, that even our relatively very little ignorance about cosmological parameters can induce an O(l) uncertainty in some 
results, since, e.g., terms like Qcdm and <rg enters to the fourth power in C'i. 
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FIG. 2: Source counts. Left: Differential number counts for AGN (dashed), star-forming galaxies (dotted), and the benchmark 
B2 annihilating DM model (solid) described in the text. Data are from Ref. [§3] and are multiplied by a factor c (given in the 
figure inset) for clarity. Central: Differential number counts at 1.4 GHz for all benchmark DM models described in Tabled] 
Right: Redshift distribution of a bright (S > 3 mjy, lower) and a faint (0.1 /iJy < S < 1 /iJy, upper) samples for astrophysical 
sources (only dominant population is shown) and the benchmark DM models B2. 



(1 + z) in the CMB target photon energy [56j . For 7r°-decay emission, there is an extra decrease with redshift which 
depends on the WIMP annihilation/decay spectrum, since the emission is given by dN 1 / dE((l + z)E). Accordingly, 
the contribution at high redshifts becomes less significant with respect to the IC on CMB case. 

An important conclusion that can be drawn from Fig. [TJ; is that the synchrotron total intensity is provided by 
emissions at low redshifts, lower than in the case of the IC case. This is easy to understand and is due to the rapid 
increase in the energy loss associated to IC scattering on CMB (scaling as (1 + z) 4 ), while a significant increase with 
z of the magnetic field is not expected. 

In Fig. [2]numbcr counts are shown. The annihilating WIMP model B2, which leads to a total intensity comparable 
to astrophysical contributions (as shown in Fig. [1]), is compared to AGN and SFG counts in Fig. [2^,. AGNs largely 
dominate at strong brightness, while SFGs take over below mJy-nuxes. The DM contribution becomes dominant in 
the sub-/xJy regime. The corresponding decaying case has a steeper spectrum (see Fig. [2b), and so, although at large 
brightness it is relatively more important than the annihilating case (but still subdominant), it takes over only well 
below the /xJy level. This is because of the p 2 scaling in the annihilating case (combined with the growing of the 
concentration parameter as the halo mass decreases) which makes the smaller and fainter structures relatively more 
important than large and bright halos. 

Fig. [2fc shows again (along the same line of Fig. [T]) that most of the WIMP-induced radio emission comes from 
sources at very low redshift. For large brightness, AGN dominates and the DM contribution is subdominant basically 
at all redshift. On the other hand, in the sub-/xJy regime, DM sources are more numerous than SFG up to z ~ 1. 
The precise shape of SFG counts is not completely understood yet, and curves in Fig. [2] represent one possible model. 
However, the overall normalization and the fact that SFGs peak at z > 1 are rather robust predictions [oTI ]. It means 
that we can firmly conclude that a WIMP source population which provides a substantial contribution to the total 
intensity has to dominate in the source count at low brightness (S < /xJy) and low redshift (z < 1). 

In Fig. [21 we show the angular correlation of radio sources. Dimensionless Ce and w are normalized to, respectively, 
the mean intensity (I) and number counts N corresponding to the astrophysical contributions. For what concerns 
the power spectrum, the one-halo term of Eq. ([5]) (which is basically a Poisson noise term, Ce ~ const, up to very 
large multipoles) dominates at large brightness (namely, at brightness currently covered by data), where the DM 
contribution is again much smaller than the AGN one. In the sub-^tJy regime, the two-halo term is instead the most 
important, and the DM population starts dominating, in particular at low multipoles, i.e., large scales. Note also 
that the DM power spectrum is completely different from a flat Ce in both Figs. [3^ and[3}3 (except when taking only 
brightest halos), telling us that it is mostly provided by the two-halo term. This is again due to the fact that the 
DM source peaks at very low redshift, which explains also why the DM tends to have more power on large scales. For 



data points corresponding to observations which integrate over all brightness, like the Haslam et al. map [57[ at 408 
MHz and the SPT survey [3{| at 150 GHz, the Ce are mostly given by the brightest objects and the DM contribution 
is only marginally relevant at I < 100 (where however any estimate of the extragalactic spectrum is undermined by 
uncertainties in the Galactic anisotropics). The decaying case is again less (more) favorable for faint (bright) fluxes 
with respect to the annihilating DM. Note also that it has relatively more power on large scales, stemming from the 
p 2 dependence of the signal in the annihilating case which enhances the relevance of small scales. A similar discussion 
applies also for the angular correlation function plotted in Fig. [3};, with slightly more favorable conclusion for WIMPs. 




FIG. 3: Angular correlation. Left: Angular power spectrum at 1.4 GHz of a bright (S > 10 mjy, lower) and a faint (S < 
1 /xJy, upper) samples for astrophysical sources (only the dominant population is shown) and the benchmark B2 annihilating DM 
model. Data points are derived from NVSS in Ref. [3lJ. Central: Angular power spectrum computed including contributions 
from sources at all bri ght ness at 408 MHz (upper, observational band derived from Ref. [13], see Refs. [13, H3]) and 150 GHz 
(lower, data from Ref. [35[) for astrophysical sources, the benchmark DM models B2 (at 408 MHz), and a WIMP model (at 150 
GHz) with M x = 1 TeV and (a a v) = 3-10 -23 cm 3 s _1 (annihilating) or r = 10 26 s (decaying) (i.e., fitting the PAMELA positron 
excess [5S|). Right: Angular correlation function at 1.4 GHz of a bright (S > 3 mjy, lower) and a faint (0.1 /iJy < S < lttJy, 
upper) samples for astrophysical (only dominant population is shown) and DM benchmark models B2. Data are from [2q . 

A big caveat in our analysis is that we have been treating the source population given by DM induced signals as 
a separate population with respect to AGN and SFG. On the other hand, any astrophysical source is embedded in a 
(typically much fainter) DM halo. This means that if one subtracts bright sources to isolate, say, sub-/xJy sources, he 
might be subtracting also emissions from the corresponding DM halos, so possibly a significant fraction of the total 
DM contribution. Therefore plots in Figs. [5] and [3] are theoretical predictions, but they could be difficult to be fully 
tested through observations, since the latter will be probably biased towards DM sources hosting a faint baryonic 
counterpart. A possible way to circumvent this issue is to compute the total emission for each structure in the universe 
(namely to add the DM and astrophysical signals arising in the same source) , and then comparing the results between 
the cases with/without DM (and applying different brightness cuts). This would require to associate to any DM 
halo of mass M and redshift z all the different possible astrophysical source population and spatial distribution in a 
statistical way. This is theoretically very challenging and deserves a dedicated analysis. 

B. Discussion on assumptions 

1. Dark matter clustering 

As described in Sec. Ill Bt we consider a 'semi-analytical' procedure where the contribution from halos of all masses 
is integrated by using analytical functional forms, which are tuned to reproduce results from numerical N-body 
simulations. The most important ingredients are the halo mass function, concentration, and density profile, the 
minimum halo mass, and the subhalo abundance. 

The halo mass function is taken from Ref. [o2l ]. Considering more recent and accurate results from numerical 
simulations, the agreement is within 10-20% (see, e.g., Fig. 7 in Ref. [Hj]). We choose the function from Ref. [H| 
since it is well behaved at all masses, while results from, e.g., Ref. [o9| would need some arbitrary extrapolation at 
small masses. 

Wc consider Refs. [6(| and [53[ for most recent results regardin g th e concentration-mass relation in the CDM 
scenario. The main difference (and most interesting feature) of Ref. [6C| with respect to Ref. [53| is that they find a 
flattening and upturn with increasing mass for halo concentrations at high redshifts. However, this is not particularly 
relevant for our purposes since the synchrotron signal is mostly generated at low redshift and in low mass halos, where 
they are in good agreement. Indeed in terms of radio emission, the two models are basically coincident as shown 
in Fig. [4] For convenience we choose the halo concentration in Ref. [53| as our benchmark parametrization since it 
is defined in terms of c V i r — M V i r (we take into account a scatter in this relation by assuming that c V i r follows a 
log-normal distribution with a width of 20% [HI). Taking an observationally-driven point of view, one could use 
the relation between the halo core radius ro and central density po found in Ref. [6l| for cored DM profiles, and 
express it in terms of c V i r (M V i r ). We found that the relation tqPo — const ~ 2 [6l[ implies halo concentrations which 
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FIG. 4: Clustering. Left: Total intensity for the benchmark annihilating model Bl considering three different halo profile, 
NFW (solid), Burkert (dashed), and Moore (dashed-dotted) , two different minimum halo mass, M cut = 1O 6 M0 (lower) and 
Mcut = 10~ 6 Mq (upper), and two different models for the concentration-mass relation, from [gO] (red) and (blue). Central: 
Number counts at 1.4 GHz for the benchmark annihilating model Bl considering an NFW profile with M cut = lO 6 A/0 (solid) 
and M cut = 1O _6 M0 (dotted), together with Burkert (dashed), and Moore (dashed-dotted) profiles (with M cu t = 1O 6 M0). 
Right: Angular power spectrum for the same models of central panel. The one (blue) and two (red) halo terms are shown 
separately. 



are similar to the results of simulations for M = 10 6 — 10 8 M Q , while significantly lower concentrations are obtained 
for larger masses. This relation stems from measurements of nearby galaxies and so we do not attempt to use it in 
our computations which involve both larger and smaller structures and halos at z > 0. However, even extrapolating 
somehow the relation in those ranges, we do not expect dramatically different predictions, since the DM signal is 
mainly given by the contribution of the smallest halos. 

A standard reference for the concentration-mass relation adopted in estimates of extragalactic signals has been 
Ref. [62| . Now it is found to be not consistent with most recent simulations [H, [6(| , mainly because underpredicts the 
concentrations for more massive halos. In our computation, however, it leads to, at most, 0(1) difference (depending 
on the choichc of K and F parameters in the model), so it does not represent a large source of difference with respect 
to past works. 

The main uncertainty comes from the extrapolation at low masses, where unfortunately we don't have much 
simulation data. In this case, different models can lead to significantly different results. To bracket such uncertainty 
we consider a naive extrapolation of the model of Ref. (53| down to 1O _6 M0 and a case where there is no contribution 
below 10 6 Mq. As shown in Fig. 3^, for the intensity there is a factor of 50 between this two pictures. Such uncertainty 
is even larger than the one associated to the halo density profile, whose effect is shown by plotting three different cases, 
NFW, Burkert and Moore. Only the latter, which is a quite extreme choice, provides a significant boost. Note that only 
the NFW profile is strictly consistent with simulations [H, [6(3, [!H adopted for halo concentration. When considering 
Burkert and Moore profiles, we arc somehow implicitly assuming that they arise from an evolution of an NFW profile 
induced by 'non-cosmologicaP effects associated to barionic physics. It's interesting to note that, for number counts, 
the choice of M cut makes, in practice, little difference. Indeed, very low mass halos have very low luminosities and 
therefore significantly increase the counts only in a range very far from experimental sensitivities. At larger luminosities 
the small difference simply comes from the normalization of dn/dM which is given by J M t dM dn/dM M = Vt m p c . 
Since we stop at a given halo mass M cut , the counts stop at brighter luminosities for steeper profiles, and the faint 
end of number counts in the Moore case drops at larger fluxes with respect to NFW and Burkert cases. A possible 
argument in favor of M cut > 10 6 M Q relies on the need of a reasonable amount of barionic matter in order to generate 
a magnetic field significantly stronger than the cosmological one (which is ~ nG). The fact that smallest DM halos 
would not (presumably) host a large amount of luminous mass component implies that synchrotron radiation would 
be accordingly much smaller than in more 'standard' halos. This will be discussed in Sections IIII B 31 and IIII B 4l 

Fig. 3J; shows the impact of halo density profile and minimum halo mass on the angular power spectrum. As for 
number counts, the choice of M cut makes little difference for the dimension- full Cp (except for mildly enhancing the 
relative importance of the 2h term with respect to lh), since the correlation of faintest halos becomes important only 
if focusing on emissions at very low brightness. Dimcnsionless Ci are strongly rescaled due to the much higher factor 
(I) (reflecting what found in Fig. 0k) in the normalization. For the 1-halo term, the Burkert and Moore profiles 
show, respectively, less and more power at small scales (I > 10 3 ) with respect to the NFW case, as expected. The 
2-halo term is instead mostly set by the linear power spectrum of Eq. |(7J) (which is obviously the same for all the 
three cases). Note also, as already mentioned above, that the 2-halo term dominates up to large multipolcs. Fig. [4]: 
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FIG. 5: Substructures. Left: Total intensity for the benchmark annihilating model Bl with (dashed/dotted) and without 
(solid) substructure contribution. The substructure model is such that 10% of total mass is in substructures with subhalo 
mass function dn B /dM s oc M~ 19 . The minimum subhalo mass is M a cut = 10 6 M Q (dashed) or M c s „ t = 10" 6 M Q (dotted). 
Central: Number counts at 1.4 GHz for the same models of left panel. Blue color stands for the case in which substructures 
are considered as single sources, while red denotes the picture in which their contribution is only a boost of the signal from the 
parent halo. Right: Angular power spectra of host-halos and sub-halos for the same two models of left panel. In the model 
with M" ut — 1O _6 M0, both biased and antibiased radial distributions of subhalos are shown. 

is obtained by integrating basically over all brightnesses (S < 1 mJy). When cutting contributions from high-flux 
sources, the 1-halo term decreases much more rapidly than the 2-halo term, which therefore can dominate the whole 
spectrum up to very large multipoles. 

Summarizing, the total uncertainty from DM clustering properties (without introducing substructures) in the 
intensity, number counts, and angular power spectrum of annihilating DM reaches about two orders of magnitude. 
Note from Eqs. (fl~8|) and ([20]) that, under the approximation of spatially constant magnetic field and diffusion which 
do not vary with the halo mass, the results in the decaying case arc completely independent from all the uncertainties 
discussed in this Section. 

To understand the possible contribution from substructures, we consider a model such that 10% of total mass is in 
substructures with subhalo mass function dn s /dM s oc M~ a and a = 1.9 which is in agreement with simulations for 
Milky Way size halos [HI, [63[ (see also similar treatments in, e.g., Ref. 0, [13]), while all other clustering details are 
kept as in the benchmark case of the previous Section. The minimum subhalo mass is taken to be either M£ ut = 10 6 Mq 
or M^ ut = 10~ 6 M Q (with the minimum mass for the parent halo fixed to M cut = 1O 6 M ). Those two models leads 
to a boost in the total intensity of a factor of 1.4 and 7, respectively, see Fig. [5^,. 

For what concerns number counts, we have to distinguish the case in which a subhalo is identified as a single source 
to the case such that the subhalo only contributes to boost the signal of the parent halo. In the first case, we say the 
subhalo is resolved, while in the second it is unresolved. To bracket uncertainties related to the possibility of resolving 
substructures, we consider two extreme cases such that subhalos are either all resolved or all unresolved. The first 
possibility leads to a flattening in the number count spectrum, since the number of sources at lower brightness is 
increased. On the contrary, the second possibility steepen the spectrum, since the boost from subhalos pushes parent 
halos to larger brightnesses. 

Fig-Efc shows the angular power spectra of the host-halos and subhalos in the scenario where both contributions are 
included 3 . The host-halo case is simply rescaled by a factor (1 — /) 4 (with / = 0.1) with respect to the plot of Fig.dk. 
With the formalism provided by Eqs. ([5]), ©, and (JTUJ), the subhalo contribution is computed in a similar fashion 
as for the main halo case, but with each main structure now "weighted" by different luminosity and spatial profile 
that follow from the (averaged) substructure distribution. The emission profile is now shallower than p 2 , and the 
emission from largest halos is enhanced since they host largest number of subhalos. The latter property implies that 
the relative significance of the one-halo term with respect to the two-halo term increases if smallest substructures are 
not introduced (case with M^ ut = 10 6 M Q ). The first property leads instead to a flattening of the power at very large 
multipoles in the one-halo term with respect to the host-halo case. The effect is visible for subhalo density going as 
oc p (biased case), and even more pronounced in the antibiased case [43| . where the radial distribution of subhalos is 



For simplicity, we only consider the autocorrelation terms of host-halo and subhalo, disregarding the cross-correlation power spectrum 
(which is important only if the two contributions are of the same order). 
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FIG. 6: Microscopic properties. Annihilating cases with a a v = 3 ■ 10~ 26 cm 3 s -1 (left), and decaying cases with r = 10 28 s 
(right). 

much natter than that of the dark matter in the host halo. 



2. Dark matter microscopic properties 

In Fig. El we show how the synchrotron signal depends on the DM microscopic properties. We consider few different 
masses, and two rather different final states of annihilations/decays, b — b and /i + — fx~, the first inducing a e~/e + 
spectrum much softer than the latter (and in turn a synchrotron spectrum with analogous properties, see the plot). 

The dependence on annihilation/decaying rate is straightforward since it is just a normalization parameter 4 . 

In the monochromatic approximation of synchrotron radiation, the energy of an electron emitting at frequency 
v is given by E ~ 15\Jvghz / B^G GeV (where vqhz is the frequency in GHz and B^g is the magnetic field in 
pG). Assuming a magnetic field in the reasonable range of few microGauss, we can get an estimate of the energy of 
electrons/positrons involved in the emission at a certain frequency. With this simple relation we can understand how 
to relate the e + — e~ spectrum of annihilation/decay with the emission. The threshold in Fig. [5] is given by the fact 
that no electrons can be produced with energy above the mass M x (in the annihilating case) or above M x /2 (in the 
decaying case), since WIMPs are non-relativistic. Then, softer is the e + — e~ spectrum of annihilation/decay, smaller 
is the energy at which the emission starts dropping. 

The normalization of the annihilation signal scales as the square of the number density, so as M~ 2 , while in the 
decaying case as M^ 1 . Therefore the spread in the decaying case is much smaller, and it is more predictive than the 
annihilating case. 

For the same mass and annihilation/decaying rate, WIMPs annihilating into p + — pT induce a signal which is 
brighter (fainter) at high (low) frequencies with respect to the b — b case. We have already stressed that the bulk of 
the DM induced emission is generated at very low redshift. Therefore there is no particular mixing among different 
redshift slices and the impact of features of the annihilation/decay spectrum on source counts and anisotropies is 
similar to what shown for the total intensity. 

The comparison between annihilating and decaying DM has been performed in Figs. Od, [3b, and[3]:. The p 2 scaling 
in the former case enhances the power on small scales and the number of faint sources with respect to the latter 
case, which scales with p. These facts makes the emission in the annihilating case more distinguishable from other 
astrophysical radio sources. 

3. Magnetic field 

First of all, we investigate the impact of the strength of the magnetic field, but still assuming B to be constant 
in time and space. The magnetic field affects both energy losses and synchrotron power. As mentioned above, we 
consider energy losses due to synchrotron and inverse Compton scattering on CMB. The case of IC scattering on CMB 



4 It could actually affect small-scale clustering, but only for very dense regions and large rate which are not considered here. 
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FIG. 7: Magnetic field. Left: Total intensity induced by the benchmark annihilating model B2 considering four different 
assumptions for the magnetic field strength. Central: Number counts at 1.4 GHz for the same models of left panel. Right: 
Angular power spectra of DM sources with brightness 0.1 /zJy < 5* < 1/^Jy at 1.4 GHz for the same models of left and central 
panels. One-halo and two-halo terms of Ct are shown in blue and red, respectively. 



is equivalent to an energy loss due to synchrotron radiation on a magnetic field of So — 3 /zG (for low energy electrons, 
i.e., neglecting Klein-Nishina corrections). Therefore for magnetic fields much larger than Bq, the dominant energy 
loss is synchrotron. In this case, the total flux is only marginally affected by the magnetic field strength. Indeed 
a decrease in B leads to a lower synchrotron power which is however compensated by the larger cooling time, and 
vice-versa. This is true provided that the confinement time is large and that the energy of electrons needed for the 
emission at a certain frequency is far from the threshold of the spectrum, which is (a fraction of) the WIMP mass. 
We have already mentioned that, for a magnetic field of few fiG, emissions at radio frequencies are mostly generated 
by electrons with energy around 1-10 GeV. For magnetic fields smaller than 3 /zG, the scaling is instead roughly 
proportional to B 2 , since only the synchrotron power is affected, while energy losses do not significantly vary. The 
case with B = 1 /iG is plotted in Fig. [TJt which shows that, with respect to the case with B = 10/iG, there is a factor 
of 10 in the normalization, and the threshold occurs for lower frequencies. 

To understand what could be the effects of a more realistic treatment of magnetic field, we consider a simple but 
instructive analytical form for B: 

B(M '"> = M£;)° (1+ * r M-i)' <22) 

with A/ max = 10 18 M Q (the exact chosen value doesn't really matter) being the largest halo mass considered. With 
this simple expression, we keep spherical symmetry and the crucial parameter to set the strength for a given source 
is again the mass M, so Eqs. B1fT5lin Section II are still valid and easily solvable (and will be exploited in the rest of 
the paper, no longer resorting to the simplified case of Eqs. [TT1I2T))) . On the other hand, Eq. (|22p allows us to study 
a magnetic field which increases with the size of the source (a > following observations of larger B in clusters with 
respect to galaxies), which has a time evolution (typically b > 0), and which decreases in the outer part of the source 
far from the luminous region (c < 1). 

The case considered so far is simply recovered by (a,b,c) — (0, 0,oo). As shown in Figs. [TJ: and the bulk of 
the DM signal is emitted at very low rcdshift. Therefore b ^ leads to negligible differences, unless considering 
extreme evolution: the z-dependence of the magnetic field can therefore be safely neglected. The scalings with object 
mass and spatial extension of B in the object can be instead more interesting and in Fig. [7] we show few examples, 
considering a = 0.1 (to go from B ~ 10 /iG in large systems to B ~ 1 /j,G in small galaxies) and c = 1/50 (e.g., in 
our Galaxy the virial radius is R v ~ 250 kpc and the magnetic region is few kpc in the vertical direction). The case 
with a > leads to a suppression in the emission from small objects, which implies a depiction in the number counts 
at small brightness (Fig. [TJ)) and an overall depletion in the total intensity (Fig. [7^). The one-halo term of angular 
power spectrum increases its relative importance with respect to the two-halo term, and the crossing point between 
the two curves in Fig. [7J; occurs at smaller I compared to the a = case. The model with c / 0, leads instead to a 
suppression in the emission from large objects. Indeed, since the concentration parameter c V i r decreases as the halo 
mass increases and thus the DM halo profile drops more quickly in smaller objects, the impact of cutting the emission 
at a certain radius <C R v is larger on massive systems. The effect on total intensity and angular correlation is however 
modest, with a very mild increase of the relative importance for the two-halo term. 
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4- Electron/positron transport in extra-galactic sources 



So far we have been considering a picture such that the e~/e + produced in WIMP annihilations/decays radiate at 
the same time and place of injection, and that all the power injected by DM is radiated. In other words, we have been 
using the solution to the transport equation for e~ /e + given by Eq. 1161 The time-scale associated to energy losses is 
taken to be: 

3.7-10 16 s GeV 
Tl ~ (l + z) 4 + 0.095 Bl G E ■ ( 3) 

We included energy losses from synchrotron and inverse Compton on CMB (here reported neglecting Klein-Nishina 
corrections for IC, which is a good approximation up to 10 TeV for scattering with CMB photons). Since we 
considered magnetic fields of 1-10 /xG and a radio emission which is associated to electrons with energy 1-10 GeV, 
we typically have 77 < 10 16 s. Three are the conditions that need to be satisfied such that Eq. (|16I23[) provides a 
good approximation of the real picture, namely 77 has to be, in each structure, significantly smaller than timescales 
associated to: 1) evolution of the structure, 2) diffusion of charged particles and 3) confinement of charged particles. 

Condition 1) is the one that allows us to adopt a time-independent solution. Since 77 is always much smaller than 
the Hubble time (as already noticed in Rcf. [56|). we can assume that e~/e + are injected at the same redshift at 
which they emit and neglect (in this respect) the cosmological evolution. Also, since the emission is concentrated at 
very low redshift (and span a very limited redshift range), we can assume e~/e + to have reached equilibrium and 
consider a steady-state solution. 

Concerning condition 2), we first notice that so far the spatial distribution of e~/e + in the source is assumed 
to follow from the DM spatial profile. This can be considered as a good approximation only if the diffusion scale 
cIl ~ a/4 D ti (where D is the diffusion coefficient D = Do (E /GeV) s ) is smaller than (comparable to) the scale at 
which the DM profile shows significant variations. If we consider the value of Do = 3 • 10 28 cm 2 /s as inferred in the 
Milky Way, we have larger than few kpc. Therefore, in particular for cuspy profiles, the spatial diffusion has an 
impact on the computed signal. It tends to smooth out the cusp of the injection distribution into shallower profiles. 
However, provided that the particles do not efficiently escape the object (see discussion below) and the magnetic field 
does not dramatically vary within the diffusive region, the brightness of single sources is not affected (and in turn also 
the total extragalactic intensity). Only the angular correlation can be modified at very small scales, in a way that 
can be understood moving from the Moore or NFW case to the Burkcrt curve in Fig. 0J:. Indeed, effectively, a cored 
profile mimics the spatial distribution of e~/e + injected with cuspy profile and then undergoing spatial diffusion. 

If condition 3) is not verified the source brightness can be depleted. Indeed it means that e~/e+ can escape the 
object before radiating at the frequency of interest. A full description of spatial diffusion in all the extragalactic sources 
is obviously beyond the goal of this paper. To estimate the impact of a finite escape length, we multiply the e _ /e + 
spectrum dN e /dE by a factor exp(— T coo i/T con f). The confinement time r COTt / can be estimated in a simple leaky box 
model as T con f(E) ~ L 2 /(2 D(E)), where L is the size of the diffusive region. The cooling time is instead t coo i(E) ~ 
Ti(E)(E/E p — 1), where E p is the electron energy corresponding to the peak in the synchrotron emission at the 
frequency of interest. Under such approximations, the confinement time in the Milky Way is r co „/ ~ 10 15 (-E/GeV) _l5 
s (having taken L to be few kpc). Therefore even a relatively large structure as our Galaxy can suffer of a depletion 
in the emission at low frequency (since the latter involves electrons with lower energies which means larger cooling 
times and diffusion lengths) due to the escape of a fraction of the injected e~/e + . Note, however, that we have 
been neglecting the energy loss due to IC on starlight which, in high-luminosity structures, can become relevant and 
reduce the cooling time. The effect of electron escape becomes more and more important for smaller objects, and 
the possible large enhancement provided by DM halos with masses down to 10 _6 M© (see Fig. 2^) can be partially 
erased by the efficient escape. This is shown in Fig. [5Jt, where we compare the total intensity of the case neglecting 
electron escape to the case where it is accounted for in the described simplified picture (for all the plots in Fig. [FJ 
we choose L = R v /50 and D = 3 • 10 28 (i5/GeV) 1//3 cm 2 /s). At frequency above 10 GHz, the effect is negligible since 
the emission mainly comes from electrons with energy larger than 10 GeV which loose energy very effectively and 
radiate nearly instantaneously after injection. Fig. \Ejp shows that number counts associated to small objects at low 
brightness arc depleted (up to extremely low brightness where the counts stem from sources at high-rcdshift having 
much shorter diffusion lengths due to the much larger IC cooling). In Fig. [8J we plot how the redshift distributions 
of large (10 6 M Q < M < 10 18 M Q ) and small (10" 6 M Q < M < 10 6 M Q ) sources are affected. In the latter case, the 
electron escape largely depletes the emission and also narrows down the contribution to lower redshift. Larger objects 
are too bright to emit at brightness 0.1 ^Jy < S < 1 /xJy for redshift z < 0.01 if the escape is neglected. When the 
latter is included, however, the contribution from sources in the lighter end of the mass range is shifted towards lower 
redshift. Most massive halos are dominant at large redshift and their emission is unaffected. 

The effect of electron escape on angular correlation at low brightness is to mildly enhance the relative power of 
large versus small scales. This is due to the depletion of the contribution from small sources and to the shift towards 
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FIG. 8: Confinement. Left: Comparison of total intensity with/without possible escape of e~/e + after injection from DM 
annihilation (in the benchmark annihilating model B2). The two cases with M c = 10 6 M© and M c — 10 -6 Mq show that 
small structures and low frequencies can be significantly affected. Central: Source counts at 325 MHz with/without possible 
escape of e~/e + and considering M c = 10 6 Mq and M c = 1O~ 6 M0 (again in the benchmark annihilating model B2). Right: 
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C. Constraints on WIMP parameter space 



Having discussed the impact of the main uncertainties on the WIMP induced signals, we can now proceed to derive 
constraints on the WIMP parameter space. They are obtained through the comparison with the data presented in 
Section 2. We conservatively choose not to add any other astrophysical contribution to the theoretical flux. 

Considering the benchmark scenario depicted in Sec. 1111 A[ we compare bounds from total intensity, source counts, 
and angular correlation in Fig. [9^. The exclusion curves scales roughly as for M x above few tens of GeV, while 
for lighter WIMPs the induced spectra can become significantly softer than data, which leads to a flattening of the 
curves. Note that the total-intensity bound is subdominant. However, as already mentioned, we have been exploiting 
extragalactic estimates from the ARCADE collaboration which are 5-6 times larger than the expected emission. 
Whether this excess will be ascribed to instrumental or galactic contaminations rather than to extragalactic sources, 
the bound should be accordingly rcscaled. Moreover, for different scenarios with an enhanced total emission due to 
the inclusion of low brightness objects (e.g., reducing M cut or introducing (resolved) substructures), the intensity 
constraint becomes more effective than number counts since for the latter there are no data at low brightness. 

As already stressed above, the experimental data currently available make constraints from anisotropics and number 
counts significantly weaker with respect to what one could in principle achieve with deeper future survey (while the 
total-intensity bound is not expected to dramatically improve). 

On the other hand, and again as discussed above, since DM halos can host much brighter baryonic components, 
the induced emission of the former can be obscrvationally hidden by the latter for what concerns measurements of 
number counts and two-point angular correlation function. We leave the discussion on how to separate the DM 
contribution for a further work. Here, to derive conservative constraints wc consider only total intensity and angular 
power spectra CV derived from observations with single dish telescopes and CMB satellites (which get contributions 
from sources at all brightnesses). The lowest multipole considered at low (high) frequency is £ = 50 (£ = 3000) because 
of contaminations from galact ic ( CMB) anisotropics which make any estimates of extragalactic source contribution 
on larger scales too uncertain |32rf34| . 

For our final constraints on WIMP parameter space, we define two sample cases bracketing the effect of quantities 
which affect the extragalactic DM radio flux but not directly related to DM microscopic properties (namely: clustering, 
substructures, magnetic field, and spatial diffusion). The first scenario (called MIN) is such that the flux is minimized, 
while the second (MAX) corresponds to the most optimistic assumptions. In both models, we adopt the halo mass 
function from Ref. [52| , the concentration of halos from Ref. [53[ , and the DM distribution in halos follows a NFW 
profile ||. 

In the MIN case, we set M cut = 10 6 M Q , no contributions from substructures, magnetic field as in Eq. (|2"2")l (with 
Bq = 10 /xG and (a, b, c) = (0.1, 0, 1/50), and electron escape as described in Sec. lIII B~4l For the MAX case, we choose 
Mcut — 10 _6 M Q , B = 10 /uG, and e + /e~ are assumed to radiate at the same place where they are injected (while the 
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FIG. 9: Constraints. Left: Bounds from total intensity (solid), source counts (dashed), and angular correlation (dotted) 
for WIMPs with dominant annihilation channel fi + — [iT . Clustering and other astrophysical assumptions are as described in 
Sec. IIII Al The reference 'thermal' annihilation rate a a v = 3 • 10 -26 cm 3 s -1 is shown for comparison. Central: Constraints on 
the parameter space of WIMP annihilating models in the MIN (solid) and MAX (dashed) scenarios (see text for details) . Final 
states of annihilations are b — b (blue) and fi + — [f (orange). Right: Same of central panel, but for decaying DM. 



inclusion of substructures makes little difference unless considering extreme models for the subhalo mass function). 

We choose again the two benchmark final states of annihilations (decays), 6—6 (inducing a softer e + — e~ spectrum) 
and fi + — /j," (inducing a harder e + — e~ spectrum), and show exclusion plot in terms of M x vs (a a v) (annihilating 
case) and M x vs r (decaying case) in Figs, and [9k- 

Although the impact of various single assumptions is quite different in annihilating models with respect to decaying 
scenarios, there is roughly three orders of magnitude between the MIN and MAX constraints in both cases. Only for 
low mass WIMPs and fairly optimistic assumptions we can close on the benchmark 'thermal' annihilation rate (a a v) = 
3 • 10 -26 cm 3 s _1 . In the decaying case, the inclusion of angular power spectrum data leads to significant improvements 
with respect to constraints derived from total intensity, especially for large masses (while in the annihilating case the 
improvement is negligible). This can be understood also from Fig. [2b- 



IV. CONCLUSION 



In this paper we presented an extensive analysis on the cosmological radio emission associated to WIMP dark matter. 
This emission is given by synchrotron radiation involving electrons and positrons injected by DM annihilations or 
decays in extragalactic structures. WIMP models with large branching ratios into light leptons lead to large e + /e~ 
yields with suppressed 7r° production and in turn no bright 7-ray emission from 7r°-decay. For those models the ratio 
between the DM signal and the background astrophysical emission produced by "standard" sources is significantly 
higher at radio frequencies than in the 7-ray band. 

For what concerns non-thermal emissions, radio telescopes have probably greater capabilities (i.e., sensitivity 
and angular resolution) as compared to telescopes at other wavelengths. This allows to study few different and 
complementary observables. We considered the total intensity, the angular correlation of source distribution and 
anisotropics, and the source number counts. We found that the emission induced by DM is mostly provided by a 
large number of small halos at redshift well below one. In general, WIMP models which provide a non-negligible 
contribution to the total extragalactic intensity take over other astrophysical sources in the number counts at low 
brightness (S < fJ-3y). There are also interesting prospects for the angular correlation at low multipoles and, again, 
very low brightness. However, it is important to stress that such conclusions are valid under the assumption that 
the source population given by DM-induced signals can be treated as a separate population with respect to other 
astrophysical sources, as, e.g., AGN and SFG. On the other hand, any astrophysical source is embedded in a (typically 
much fainter) DM halo which makes the separation of the DM population awkward. 

We discussed the impact on the expected signals of dark matter clustering, magnetic fields, spatial diffusion, and 
dark matter microscopic properties. We found that bounds in the WIMP parameter space are uncertain by roughly 
three orders of magnitude. Under optimistic assumptions, they close on light WIMPs with thermal annihilation cross 
section, and, in general, are competitive with constraints from other indirect dark matter searches. 

The extragalactic background induced by DM can account for the ARCADE excess [17[ and the faint edge of 
differential source counts and angular correlation data, which will be provided by the next -generation radio surveys, 
can become an important probe for WIMPs. 
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